Generalized BPS magnetic monopoles 
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We show the existence of Bogomornyi-Prasad-Sommerfield (BPS) magnetic monopoles in a gener- 
aUzed Yang-Mills-Higgs model which is controlled by two positive functions, g ((^"f/)") and / {(j)^ (j)"") . 
This effective model, in principle, would describe the dynamics of the nonabelian fields in a chromo- 
electric media. We check the consistency of our generalized construction by analyzing an explicit 
case ruled by a parameter jS. We also use the well-known spherically symmetric Ansatz to attain 
the corresponding self-dual equations describing the topological solutions. The overall conclusion 
is that the new solutions behave around the canonical one, with smaller or greater characteristic 
length depending on the values of /3. 

PACS numbers: ll.lO.Kk, ll.lO.Lm, ll.lO.Nx 



I. INTRODUCTION 

In the context of classical field theories, configurations 
with nontrivial topology, generically named topological 
defects, arise as finite energy solutions to some nonlin- 
ear models. In the standard approach, such models are 
usually endowed by a symmetry breaking potential for 
the matter self-interaction, since topological defects are 
known to be formed during symmetry breaking phase 
transitions 

The simplest topological defect is the static kink, which 
arises in a (l-|-lY-dimensional theory containing a single 
real scalar field Other well-known examples of topo- 
logical structures are the vortex configurations, which 
appear in effective planar models containing a complex 
scalar field coupled to an abclian gauge field Q , and the 
magnetic monopole, which stand for the static solution 
arising from a (l-|-3)-dimcnsional theory describing the 
interaction between a real scalar triplet and non-abelian 
gauge fields [I]- 

In particular, magnetic monopoles are spherically sym- 
metric configurations coming from a static Yang-Mills 
theory endowed with a fourth-order Higgs potential @. 
In this case, these solutions exhibit no divergences and 
possess finite total energy, since they are constrained by a 
set of suitable boundary conditions. On the other hand, 
in the absence of a Higgs potential, magnetic monopoles 
arise as the minimal energy configurations of the cor- 
responding Yang-Mills model 0. In this context, such 
solutions come from a set of first-order differential equa- 
tions, and minimize the total energy of the overall theory. 

Furthermore, during the last years, some generalized 
or effective field theory models, generically named k- 
ficld theories, have been intensively investigated. The 
main difference between them and their canonical coun- 
terparts is the presence of nonstandard kinetic terms, 
which change the dynamics of the overall system in an 
exotic way. In fact, such theories have been used as ef- 
fective models mainly in Cosmology, with the so-called 
k-essence models suggesting new insights about the 
accelerated infiationary phase of the universe [§[. The 



interesting point is that the introduction of generalized 
kinetic terms has important consequences on the forma- 
tion of topological structures. Despite the possibility of 
achieving nontrivial topologically configurations even in 
the absence of the spontaneous symmetry breaking phe- 
nomenon some k-field models endowed with spon- 
taneous symmetry breaking also support topological de- 
fects In the last case, the resulting solutions can 
be studied via the comparison between them and their 
canonical counterparts, with the observation that they 
usually present slight variations on some of their main 
features 11 1. 

On the other hand, in the presence of a generalized 
dynamics, the resulting model is highly nonlinear, its so- 
lutions being quite hard to find, even in the presence of 
suitable boundary conditions. In this context, the de- 
velopment of a consistent first-order framework is quite 
useful, since new topological configurations can be found 
by solving a set of generalized BPS equations. As it hap- 
pens in the usual case, these new configurations minimize 
the energy of the overall system by saturating its lower 
bound [12| . 

Recently, some of us have performed the development 
of generalized first-order frameworks regarding several ef- 
fective field theories [l^. However, for simplicity, these 
models were constructed containing only abclian fields. 
In the present work, we go further by presenting a first- 
order theoretical framework consistent with a general- 
ized Yang-Mills-Higgs model. This model is controlled 
by two dimensionless positive functions, g {4'°'(j)°') and 
/ {(/)"■ (j)'^) , which change the dynamics of the non-abelian 
fields in a nonusual way. Nevertheless, in order to guar- 
antee the self-duality of the model, g and / arc related by 
means of a simple constraint. As in the usual case, the to- 
tal energy of the generalized model is also bounded from 
below, and the spherically symmetric self-dual solutions 
describe generalized BPS magnetic monopoles. More- 
over, in the appropriate limit, our framework leads to 
the usual one, as expected. 

In order to present our results, this work is outlined as 
follows. In Sec. |TT1 we introduce the nonstandard Yang- 
Mills-Higgs model and the corresponding first-order the- 
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orctical framework, including its generalized BPS equa- 
tions. Then, instead of recovering the usual case, we in- 
troduce an explicit generalized model, which is controlled 
by a single real parameter, /3. In Sec. IIIIl we present the 
numerical solution of the generalized BPS equations ob- 
tained by means of the relaxation technique. We depict 
the profiles of the Higgs and non-abclian gauge fields de- 
scribing the generalized BPS magnetic monopoles. We 
also comment on the main features of these solutions. 
In Sec. IIV[ we finalize by summarizing our results and 
stating our perspectives. 



II. THE THEORETICAL FRAMEWORK 

We begin introducing a generalized version of the 
Yang-Mills-Higgs model which is defined by the following 
Lagrangian density: 



C = - 



4 " 



(1) 



where all fields, coordinates and parameters are sup- 
posed to be dimensionless (this can be achieved by using 
the appropriate mass rescaling transformations) . We use 
standard conventions, including the natural units system. 
Here, 



is the usual non-abelian field strength, and 



(3) 



stands for the non-abelian covariant derivative, i]'^'^ ~ 
(H ) is the metric for the (l-|-3)-dimcnsional space- 
time, and e°'"^ is the completely antisymmetric tensor 
(with £^23 ^ -yy^g pQj^|- ^^^g^^ g {(j)" (t)") and 

f [(jf'cjf') are arbitrary functions which change the dy- 
namics of the overall model in a nonusual way. 

It is well-known that, in the usual case (i.e., given 
g = / = 1), self-dual configurations only exist in the 
absence of the Higgs potential. So, for simplicity, we 
also assume that their generalized counterparts only ex- 
ist for F (0°(/)") = 0. Furthermore, we suppose that such 
generalized configurations are described by the standard 
spherically symmetric Ansatz 



and A^ = eiai 



{W{r)-l), (4) 



where 



x°'Xa- Since we are searching for static config- 
0. The profile functions H (r) 
and W (r) are supposed to obey the usual finite energy 
boundary conditions 



urations, we choose A^ 



H{0)=0 and W (0) = 1, 



iJ(oo) = Tl and I^(oo) = 0. 



(5) 



(6) 



These conditions guarantee not only the existence of fi- 
nite energy configurations, but also the breaking of the 
S0(3) symmetry inherent to the model ([T|). 

In general, the weight functions g {(j)"' (j)"') and / ((/)°^") 
can be arbitrarily chosen, so that the Euler-Lagrange 
equations for H (r) and W (r) are quite hard to solve, 
even in the presence of the suitable boundary conditions 
([5]) and In this context, the development of a con- 
sistent self-dual theoretical framework is more than de- 
sirable, since it allows to find a set of first-order differen- 
tial equations that leads to finite energy configurations 
attained by numerical procedures. These configurations 
are genuine solutions of the overall model, once they au- 
tomatically fulfill the nonstandard Euler-Lagrange equa- 
tions coming from Lagrangian ([T]). 

From now on we focus our attention on the develop- 
ment of a BPS theoretical framework consistent with the 
nonstandard model ([T]). In order to perform it, we follow 
the usual approach, observing that the first-order equa- 
tions come from the minimization of the total energy of 
the system, given by the energy-momentum tensor zero- 
zero component. In the present case, such tensor reads 

Txp = - gF^^F^^^ - VXpC, (7) 

from which one gets the nonstandard energy density (al- 
ready written in terms of H (r) and W (r)) 



27- 



1 fdnV f HwV^ 



+^ Uu- 



(8) 



where g = g {H) and f ~ f (H) . It is clear that such 
functions must be positive in order to avoid problems 
with the energy of the model. 

Given the Ansatz (jlj, we point out that this model 
only engenders self-dual solutions when g and / arc re- 
lated to each other as 



9iH) 



1 



fiH) 



(9) 



In this case, the energy density (|8]) can be rewritten in 
the form 



2 \ dr er'^ f 



1 



/ dW 



TefHW 



(10) 



er^ dr 

whose minimization leads to the following first-order 
equations: 



dH _ l-W^ 
dr ~'~ 



dW 
dr 



er^f ' 



±efHW. 



(11) 



(12) 
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Relations (|TT|) and (|T2)) are the self-dual (BPS) equations 
of the model ([1]). After implementing these equations, the 
BPS energy density becomes 



and the total energy of the solutions is given by 

47r 



Ebps = 47r 



Sbpsdr 



(13) 



(14) 



whenever the boundary conditions ([5]) and ([6|) are ful- 
filled. 

In summary, for a given positive function f{H), the 
first-order equations (jlip and (|12p must be numeri- 
cally solved in accordance with the finite energy bound- 
ary conditions (O and ©• The self-dual solutions 
achieved in this way describe the generalized BPS mag- 
netic monopoles arising from this nonstandard Yang- 
Mills-Higgs model ([T]), with total energy given by (fT4|) 
and energy density by (|13p . However, it is worthwhile 
to point out that the generalized first-order framework 
presented in this letter only holds in the absence of the 
potential V {(j>°'(f>°'). Furthermore, g and / have to obey 
the constraint On the other hand, for a nonvanish- 
ing potential, or in the absence of the relation such 
framework does not hold anymore, and the corresponding 
BPS monopoles can not be achieved. 

The usual results are trivially recovered by setting 
/ = 1. In order to show how this generalized framework 
works, we adopt the following generalization function: 



fiH) = {H' + l) 



(15) 



where /? is some real number; here, /3 = leads us back 
to the canonical model. Now, given (fT5|) . the self-dual 
expressions PT|) and ([T^ become 



dH 
dr 



1 - 



er2 (ij2 + 1) 



= ±eHW {H^ 

dr 



1 



(16) 



(17) 



which must be solved respecting the conditions ([5]) and 

In the next Section, we solve the first order equations 
p6)) and (fTT)) by means of the relaxation technique for 
different values of /3. Then, we plot not only the numer- 
ical results for H (r) and W (r) , but also those for the 
BPS energy density p^ . and for r'^Sbps (the integrand 
of m])). We also comment on the main features of the 
new solutions. 



III. NUMERICAL RESULTS 

Now, we focus our attention on the examination of 
the profiles of the generalized BPS solutions. Thus, we 




FIG. 1: Solutions to H (r) for /3 = —2 (dot-dashed green line), 
P — —1 (dashed blue line), /3 = (usual case, solid black line), 
P — 1 (long-dashed orange line) and /3 = 2 (dotted red line). 



numerically solve the first-order equations (fT6|) and (fT7|) 
obeying the finite energy boundary conditions ([5]) and 
Here, for simplicity, we choose e = 1, and consider 
only the lower sign in dH), ([16]) and Then, 

we numerically solve the self-dual system by means of 
the relaxation technique, for different values of the real 
parameter /3. 

Note that (3 = leads us back to the usual model, 
whose self-dual solutions (already written according our 
conventions) can be attained analytically as 



Hir) = 



1 



tanh r 



W{r) 



sinh r 



(18) 



(19) 



The numerical solutions we found for H (r) and W (r) are 
depicted in Figs. 1 and 2, for (3 = —2 (dot-dashed green 
line), /3 = — 1 (dashed blue line), (3 = 1 (long-dashed 
orange line) and /? = 2 (dotted red line). The usual 
(analytical) profiles are also shown (solid black line), for 
comparison. Moreover, we also plot the corresponding 
solutions for the BPS energy density ([T^ and for r'^Sbps 
(the integrand of ([H])) in Figs. 3 and 4. 

In Figure 1, we present the numerical results regard- 
ing the profile of the function H (r) . In this case, we 
clearly see that the nonstandard solutions turn out de- 
picted around the usual counterpart, behaving in the 
same general way, but changing the width of the defect. 

We point out that different solutions engender different 
characteristic lengths. In particular, the solutions aris- 
ing for /3 < reach the asymptotic condition H {go) = 1 
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FIG. 2: Solutions to W (r). Convention as in FIG. 1. FIG. 4: Solutions to r^etpa- Convention as in FIG. 1. 




FIG. 3: Solutions to etpa- Convention as in FIG. 1. 



slower than the standard profile. In this sense, these so- 
lutions possess a greater characteristic length, so that the 
corresponding bosons mediate long-ranged interactions. 
On the other hand, the solutions related to /3 > reach 
the saturation region faster, exhibiting smaller character- 
istic lengths, which is associated with small-ranged inter- 
actions. 

In Fig. 2, we plot the solutions for W (r). Again, the 
nonstandard solutions behave as the canonical one, ex- 
hibiting variations concerned with the defect width and 



the characteristic length. Such variations are controlled 
by the real parameter /3 in the same way as before: for 
/3 < 0, the nonusual solutions engender the greater char- 
acteristic lengths, since they reach the asymptotic con- 
dition W (oo) = more slowly, being associated with 
long-ranged mediating bosons. Furthermore, for and in- 
creasing positive /3, the solutions reach the saturation re- 
gion faster, exhibiting smaller and smaller characteristic 
lengths, with the mediating bosons yielding small-ranged 
interactions. 

The numerical profiles found for the BPS energy den- 
sity Ebps (fT3)) are depicted in Fig. 3, revealing that such 
solutions are lumps centered at r = 0. Also, they vanish 
monotonically as r goes to oo (in fact, Sbps [r oo) — > 
arises in a rather natural way from the asymptotic bound- 
ary conditions ([6])). Here, the parameter /3 plays a dif- 
ferent role, since it changes not only the characteristic 
lengths of the corresponding solutions, but also their am- 
plitudes. In this sense, the solutions corresponding to 
/3 < (/3 > 0) achieve the smaller (larger) amplitudes. 

Finally, we depict the solutions for r^Sbps-, i-e., the in- 



tegrand of p4)) : see Fig. 4. In this case, the solutions are 
rings centered at r = 0. Moreover, the points of larger 
amplitudes are located at some finite distance R from 
the origin (in this sense, R stands for the "radius" of the 
ring) , such amplitudes being controlled by /? in the same 
way as before. Here, we point out the existence of an 
interesting compensatory effect: the solutions reaching 
the greater amplitudes spread over smaller distances, and 
vice-versa. As a consequence, different solutions enclose 
the same area (equal to the unity, according our conven- 
tions), and the corresponding configurations achieve the 
same total energy; see (H^ 
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IV. ENDING COMMENTS 

In this work, we have presented a nonstandard first- 
order framework consistent with a generaUzed Yang- 
Mills- Higgs model (HI, such model being controlled 
by two dimensionlcss functions, g {(j)"' (jf") and f {(jf'cjf'), 
which change the dynamics of the non-abelian fields in 
a non-usual way. In order to avoid problems with the 
energy of the system, these functions are supposed to be 
positive. We point out that the self-duality of the overall 
model only holds when g and / are related to each other 
by a simple constraint; see ([9]). Also, it is worthwhile to 
note that there is no additional constraint to be imposed 
on /. 

The generalized first-order framework was developed 
in a general way. So, in order to verify the consistency of 
our construction, we have introduced an explicit example 
controlled by a single real parameter /3; see ([T5|). Then, 
we have considered spherically symmetric self-dual con- 
figurations, the non-abelian fields being described by the 
standard static Ansatz Moreover, the profile func- 
tions H (r) and W (r) were supposed to obey the usual 
finite energy boundary conditions given by ([5]) and ©. 

The resulting first-order equations were solved numer- 



ically by means of the relaxation technique, and the self- 
dual solutions we found were plotted in Figs. 1, 2, 3 
and 4, for different values of j3. The standard analytical 
solutions (achieved for /3 = 0) were also plotted, for com- 
parison. The overall conclusion is that the nonstandard 
solutions behave in the same general way the usual one 
docs, the main difference being slight variations on the 
amplitudes and on the characteristic lengths of the new 
solutions. In addition, we have identified the way such 
variations are controlled by the real parameter /?. 

Recently, some of us have performed detailed investi- 
gations addressing generalized self-dual frameworks for 
abclian models, attaining their respective first-order so- 
lutions [l^. This way, the present letter is an extension 
of those works to the non-abelian context. As for future 
investigations, interesting issues including the supersym- 
metric extension of the non-abelian model ([1]), and the 
search for its topological structures, are now under con- 
sideration. 

The authors thank Caroline dos Santos for useful com- 
ments about this work. We arc also grateful to CAPES, 
CNPq and FAFEMA (Brazil) for partial financial sup- 
port. 
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